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Evaluation  of  a  Portable  Electromagnetic  Induction 
Instrument  for  Measuring  Sea  Ice  Thickness 


AUSTIN  KOVACS  AND  REXFORD  M.  MOREY 


INTRODUCTION 

Sea  ice  is  a  multi-component  medium.  Its  constituent 
parts  are  fresh  ice,  liquid  brine  inclusions,  gas  pockets 
and,  depending  upon  eutectic  factors,  solid  salt  crystals. 
The  volume  of  fresh  ice  is  by  far  the  larger  fraction, 
typically  in  excess  of  95%.  Sea  ice  is  classified  by  age 
(first-year,  second-year  and  multi-year)  and  by  mor¬ 
phology.  Variations  in  growth,  melt  and  deformation 
processes  result  in  ice  formations  of  complex  shape, 
structure,  and  brine  and  gas  contents.  In  particular,  the 
complex  structure  and  liquid  inclusion  variations  have 
greatly  limited  our  ability  to  remotely  measure  sea  ice 
thickness.  This  is  especially  true  for  sounding  systems 
operating  at  VHP  frequencies  and  above.  At  these 
frequencies,  the  propagation  of  electromagnetic  energy 
in  sea  ice  suffers  high  attenuation  as  a  result  of  the 
conductive  brine  inclusions  that  increase  in  volume 
with  depth  (Kovacs  et  al.  1987a).  Kovacs  et  al.  (1987a) 
showed  that  the  conductivity  of  sea  ice  varies  with  ice 
depth,  increasing  temperature  and  brine  volume.  How¬ 
ever,  the  bulk  dc  conductivity  of  Arctic  sea  ice  seldom 
exceeds  0.05  S/m.  Only  during  the  early  part  of  the  melt 
.season,  when  the  solid  salts,  which  precipitated  out 
during  the  cold  winter  months,  redissolve  to  increase 
the  sea  ice  brine  volume,  might  the  bulk  conductivity 
reach  about  0.07  S/m. 

The  capability  of  remotely  measuring  homogeneous 
sea  ice  thickness,  to  a  high  degree  of  accuracy,  using  a 
hand-held  instrument  has  long  been  desired  by  those 
needing  to  make  quick  assessments  of  sea  ice  bearing 
''apacity  for  aircraft  and  vehicle  operations  on  the  Arctic 
Ocean  pack  ice.  Many  techniques  have  been  tried,  most 
with  limited  success  because  of  the  sea  ice  brine  content 
and  related  ice  conductivity.  Most  of  these  devices 
would  not  qualify  as  being  hand-held  and  highly  por¬ 
table  except  for  electromagnetic  induction  (PMI) 
sounding  equipment.  Since  these  EMI  systems  gener¬ 
ally  operate  in  the  VLF  frequency  band,  the  measure¬ 


ments  are  significantly  less  affected  by  the  relatively 
low  bulk  conductivity  of  the  sea  ice  than  systems 
operating  in  higher  frequency  bands. 

This  report  discusses  the  results  obtained  in  April 
1 990  using  a  Geonics  EM3 1  -D  (henceforth  called  EM3 1 ) 
terrain  conductivity  measurement  instrument  with  a 
plug-in  processor  module  that  provides  a  digital  display 
of  sea  ice  thickness.  Also  discussed  is  a  simple  method 
for  using  just  the  EM3 1  's  conductivity  measurement  to 
determine  sea  ice  thickness. 


EM31  SOUNDING  CONCEPTS 

The  Geonics  EM31  is  a  9-kg,  man-portable  instru¬ 
ment  designed  to  measure  apparent  ground  conductiv¬ 
ity  by  means  of  electromagnetic  induction  (Geonics, 
Ltd.  1984).  It  has  a  transmit  (Tx)  coil  and  a  receive  (R^) 
coil  that  function  as  magnetic  dipole  antennas.  The  coils 
are  spaced  3.66  m  apart  at  each  end  of  a  tubular  support. 
For  sea  ice  thickness  sounding,  the  coils  were  mounted 
vertically  co-planar  and  therefore  functioned  as  hori¬ 
zontal  dipole  antennas.  Sea  ice  is  relatively  resistive  and 
thus  quite  transparent  at  the  EM3rs  operating  fre¬ 
quency.  Therefore,  during  sea  ice  sounding  the  trans¬ 
mitted  (primary)  electromagnetic  field  induces  eddy 
currents  primarily  in  the  conductive  seawater.  These 
currents  in  turn  produce  a  secondary  electromagnetic 
field  that  is  sensed,  along  with  the  primary  electromag¬ 
netic  field,  by  the  receiver  coil.  The  EM3 1  is  designed 
to  measure  the  magnitude  of  the  in-phase  and  quadra¬ 
ture  components  of  the  secondary  magnetic  field.  These 
components  are  normalized  by  dividing  them  by  the 
magnitude  of  the  received  primary  electromagnetic 
field  component.  Given  that  .sea  ice  is  relatively  trans¬ 
parent  at  9.8  kHz,  the  response  measured  by  the  EM3 1 
will  be  a  strong  function  of  the  instrument  height  above 
and  the  conductivity  of  the  seawater.  Therefore,  accu¬ 
rate  measurement  of  the  electromagnetic  field  response 
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from  the  seawater  and  a  full  solution  analysis  of  the  data 
using  the  numerical  procedure  of  Anderson  (1979) 
shoi'ld  provide  a  good  estimate  of  instrument-seawater 
distance,  or  the  ice  thickness,  when  the  EM3 1  is  resting 
on  ice  of  uniform  thickness. 

The  measured  EM31  response  is  not  a  point  mea¬ 
surement  but  an  integrated  depth-volume  measurement 
w  ith  a  ciuasi-circular  footprint.  The  footprint  diameteris 
on  the  order  of  two  or  three  times  the  instrument’s 
height  above  the  conductive  seawater  surface,  depend¬ 
ing  on  Tx-Tr  coil  orientation  (Kovacs  et  al.  1987b,  Liu 
and  Becker  1990).  Therefore,  as  the  instrument  is  el¬ 
evated  above  the  seawater,  the  footprint,  about  three 
times  the  EM3 1  height  above  seawater  for  the  vertical 
co-planar  coil  orientation,  increases.  For  ice  of  rela¬ 
tively  uniform  thickness,  this  should  not  pose  a  problem, 
but  on  ice  w  ith  appreciable  undulating  bottom  relief,  for 
e.xample,  the  bottom  of  most  multi-year  sea  ice,  the 
resulting  instrument-determined  ice  thickness  will  be 
an  ''average"  one  for  an  area  around  the  instrument. 
That  is.  the  ice  thickness  measured  through  a  drill  hole 
directly  below  the  instrument  on  multi-year  ice  will  not 
likely  agree  with  the  ice  thickness  estimated  from  the 
EM3 1  measured  response. 

PREVIOUS  E.VII  STUDIES 

An  early  trial  using  two  types  of  portable  EMI 
instruments  was  conducted  by  Sinha  ( 1976).  While  he 
did  experience  equipment  calibration  problems,  he  was 
able  to  demonstrate  the  potential  of  small,  lightweight 
EMI  equipment  for  measuring  sea  ice  thickness.  This 
study  was  followed  by  those  of  Hoekstra  et  al.  (1979) 
and  Hoekstra  ( 1 980).  The  former  employed  a  Geonics 
EM31  instrument  with  an  operating  frequency  of  39.2 
kHz  and  the  latter  study  used  an  EM3 1  instrument  with 
a  now  standard  operating  frequency  of  9.8  kHz. 

Hoekstra  et  al.  (1979)  tested  the  instrument  above 
saline  ice  grown  in  an  outside  test  basin  and  then  on  sea 
ice  in  Mackenzie  Bay.  Canada.  The  test  basin  measure¬ 
ments  were  very  encouraging  in  that  good  correlation 
was  found  between  the  measured  in-phase  re.spon.se and 
the  instrument  elevation  above  the  ice.  This  led  to  the 
arctic  field  trial,  which  was  not  particularly  rewarding. 
For  example,  on  2-m-thick  sea  ice,  the  ice  thickness 
determined  from  the  instrument  measurements  varied 
up  to  4()fc.  Another  field  study  was  made  on  the  sea  ice 
in  Stefansson  Sound,  located  northofDeadhorse,  Alaska 
( Hoekstra  1980).  Only  seven  measurements  were  made, 
four  on  ice  1 .70  m  thick.  The  inferred  ice  thickness  at  the 
1 .70-m-thick  ice  sites  varied  from  about  1 .83  to  2.20m. 
While  these  limited  results  were  not  especially  good, 
they  did  indicate  that  if  the  instrument  was  properly 


calibrated  and  the  measured  respon.se  properly  evalu¬ 
ated,  homogeneous  first-year  sea  ice  thickness  should 
be  reasonably  detemnnable. 

Further  tests  using  the  EM3 1  to  estimate  sea  ice 
thickness  were  made  by  the  oil  industry.  The  results 
have  not  been  leleased.  However,  as  a  result  of  his 
knowledge  of  these  tests,  D.C.  Echert  at  Flow  Research, 
Inc.,  pursued  further  evaluation  of  the  EM3 1  for  the 
remotemeasurement  of  .sea  ice  thickness.  This  included 
both  desk  studies  and  field  trials  off  the  Alaskan  Beau¬ 
fort  Sea  coast  near  Prudhoe  Bay  (Echert  1986).  The 
results  showed  the  advantage  of  using  a  vertical  co- 
planar  (versus  horizontal  co-planar)  coil  configuration 
and  demonstrated  that  with  this  coil  arrangement  first- 
year  sea  ice  thickness  could  be  estimated,  generally 
with  a  deviation  of  less  than  15%  from  the  drill-hole- 
measured  thickness.  This  coil  orientation  was  u.sed  in  all 
subsequent  EM3 1  sea  ice  thickness  measurement  studies. 
Like  Hoekstra  (1980).  Echert  used  both  the  in-pha.se 
and  quadrature  phase  components  of  the  received 
magnetic  field  for  estimating  thickness.  These  results 
were  obtained  with  the  in-phase  of  an  EM3 1  instrument 
zero  calibrated  over  highly  resistive  pemiafrost.  The 
zero  level  of  the  quadrature  phase  is  set  by  Geonics  and 
does  not  need  recalibration.*  Echert  indicated  that  if  the 
instrument  had  been  calibrated  over  a  known  thickness 
of  sea  ice,  significantly  better  ice  thickness  estimates 
may  have  resulted. 

Because  of  these  favorable  results,  further  evalua¬ 
tion  of  EMI  sounding  ensued  (Echert  et  al.  1989).  The 
main  objective  was  to  provide  additional  internal  data 
processing  capability  to  a  standard  9.8-kH2  EM31  in¬ 
strument  that  would  enable  direct  numerical  display  of 
ice  thickness.  This  was  achieved  through  the  u.se  of  a 
look-up  table,  using  the  full  solution  multilayer  analysis 
of  Anderson  (1979),  as  provided  in  Geonics  program 
PCLOOP.  and  an  interpolation  algorithm.  This  ap¬ 
proach  assumes  that  the  in-phase  and  quadrature  com¬ 
ponents  of  the  received  magnetic  field  are  unique  to 
specific  sea  ice  thickness  and  sea  ice  and  seawater 
conductivities.  The  Flow  Research  look-up  table  was 
developed  using  10  mS/m  for  thi*  bulk  conductivity  of 
the  .sea  ice,  a  seawater  conductivity  range  from  2  to  3  S/ 
m  in  0.25-S/m  increments,  and  a  sea  ice  thickness  range 
from  0.25  to  6.0  m  in  0.25-m  increments.  The  ice 
thickness  displayed  is  an  interpolation  between  the 
tabulated  data  and  the  measured  EM3 1  response. 

Field  testing  of  the  reconfigured  EM3 1  was  done  in 
the  spring  of  1989  on  sea  ice  north  of  Prudhoe  Bay, 
Alaska  (Echert  et  al.  1989).  Sea  ice  between  0.4  and 
about  3.2  m  thick  was  measured.  For  ice  over  I  m  thick. 


^Personal  communicalion  with  J.D.  VIcNcill.  Geonics. 
[.Id..  1990, 


T 


theEM3 1  and  processormodule  (PM)  system  estimated 
ice  thickness  with'"  about  59)  of  the  drill-hole-mea¬ 
sured  thickness.  However,  the  instrument-detemiineJ 
values  became  progressively  less  accurate,  less  than  the 
direct  tape-measured  values,  with  decreasing  ice 
thickness  below  1  m.  Echert  et  al.  (1989)  suggested  that 
this  may  have  resulted  from  using  a  constant  bulk  sea  ice 
conducti  V  ity  in  the  construction  of  the  look-up  table  that 
did  not  adequately  address  the  higher  bulk  conductivity 
that  can  be  expected  for  thinner  sea  ice  (Kovacs  et  al, 
1987a).  However,  other  factors  may  have  affected  the 
results,  such  as  improper  instrument  calibration,  inap¬ 
propriate  look-up  tables,  etc. 

The  geiicrally  favorable  ice  measurements  obtained 
with  the  modified  EM3 1  -PM  instrument  by  Echert  et  al. 
( 1989)  indicated  that  this  device  may  prove  useful  for 
gathering  ice  thickness  infomiation  during  our  contin¬ 
ued  evaluation  of  airborne  electromagnetic  induction 
sounding  technology  for  tiie  measurement  of  sea  ice 
thickness  (Kovacs  et  al.  1987b,  1989).  Therefore,  the 
modified  EM3 1-PM  instrument  was  obtained,  from  G. 
White  of  Flow  Research  in  Kent,  Washington,  on  a  trip 
to  Alaska  in  April  1990  for  our  study. 

On  the  day  the  EM3 1-PM  was  pick>.d  'jp.  White  gave 
us  a  quick  review  on  how  to  calibrate  and  operate  the 
instrument.  In  addition  a  brief  operations  text  was 
provided  for  future  field  use.  At  the  time  we  were  not  in 
po.ssession  of  Hoekstra’s  orSinha's  papers  (previously 
cited),  describing  their  use  of  EMI  sounding  for  measur 
ing  sea  ice  thickness,  nor  did  we  have  the  cited  reports 
of  Echert  with  us.  We  went  into  the  field  to  use  the 
EM3 1  -PM  as  a  fully  developed  operational  instrument 
for  sea  ice  thickness  measurement  and  to  detennine  if  it 
could  provide  thicknesses  that  were  within  5%  of  the 
direct  drill-hole-measured  values  as  needed  for  our 
study. 

BEAUFORT  SEA  FIELD  TRIALS 

The  F.M31-PM  instrument  was  first  used  on  first- 
year  sea  ice  1 .6  to  1 .7  m  thick.  Here,  the  instrumen*  was 
calibrated  at  a  ite  ot  known  snow  plus  sea  ice  thickness 
as  detennined  by  a  drill  hole  measurement.  It  should  be 
pointed  out  that  at  all  our  measurement  sites  the  snow 
cover  w  as  not  removed.  Therefore,  the  measured  thick¬ 
ness  discussed  hereafter  is  that  of  the  combined  snow 
plus  sea  ice.  In  the  instrument  calibration  process,  the 
known  seawater  conductivity  was  used,  which  was  2.5 
.S/m  for  our  study  area. 

On-ice  calibration  of  the  EM3  i  -PM  in.^trument  is  the 
most  accurate  and  was  the  method  used.  This  procedure 
required  that  the  instrument  be  elevated  at  two  different 
heights  above  the  seawater.  At  each  elevation,  the 


distance  to  the  seawater,  as  determined  with  a  drill  hole 
and  tape  measurement,  was  manually  enteied,  via 
thumbwheel  dials,  into  the  PM.  The  seawater  conduc¬ 
tivity,  2.5  S/m,  was  .o  entered  into  the  PM's  memory 
via  the  thumbwheel  dials.  The  instrument  would  then  be 
operated  and  the  measureu  response  used  internally  to 
match  up  with  the  Flow  Research  look-up  table  re- 
.spop.se  values  and  the  given  instrument-to-seawater 
distance.  Through  this  process,  the  PM's  look-up  table 
was  calibrated  against  the  measured  conditions. 

After  the  instrument  was  calibrated,  it  was  taken  to 
various  locations  on  the  first-year  sea  ice  floe  where  a 
measurement  was  made  with  the  instrument  resting  on 
the  surface.  In  addition,  at  one  site  the  instrument  was 
elevated  to  heights  up  to  1  m  above  the  surface  and 
soundings  made.  At  all  measurement  sites  on  this  ice 
fioe,  the  instrument  provided  distances  to  the  ice/water 
interface  that  were  within  0. 10  m  of  the  measured 
distance. 

The  instrument  was  next  taken  to  a  nearby  area  with 
0.17  m  of  sea  ice.  With  the  instrument  resting  on  the 
surface,  a  reading  of  ice  thickness  was  made.  This 
reading  indicated  0.77  m  of  sea  ice  or  0.6  m  more  than 
existed. 

To  determine  what  the  instrutnent's  lower  ice  thick¬ 
ness  limitation  might  be.  it  was  elevated  in  increments 
above  the  surface  by  resting  it  on  cardboard  boxes.  A 
plot  of  the  EM31-PM  instrument-determined  distance 
to  the  seawater  versus  the  tape-measured  distance  is 
.shown  in  Figure  1.  This  figure  suggests  that  ice  thick¬ 
ness,  or  distance  to  .seawater,  of  less  than  about  0.7  m 
cannot  be  measured  to  within  ±10%  of  the  true  distance 
whcii  the  EM31-PM  instrument  is  resting  on  the  ice 
surface.  However,  thin  ice  should  be  adequately  mea- 
.surable  if  the  instrument  is  elevated  0.7  m  or  more  above 
the  ice  surface  and  this  distance  is  then  subtracted  from 
the  instrument-detennined  distance  to  the  seawater.  Of 
course,  this  requires  the  operator  to  judge  when  he  may 
be  on  ice  less  than  0.7  m  thick  and  must  elevate  the 
instrument.  For  relatively  accurate  EMI  ice  thickness 
sounding,  the  EM3 1  should  not  be  held  at  waist  height 
while  making  routine  soundings,  as  variations  in  the 
instrument's  coil  from  a  vertical  co-planar  orientation 
will  cau,se  measured  response  variations  that  cannot  be 
properly  analyzed  by  the  PM. 

The  results  in  Figure  1  also  show  some  scatter 
between  the  tape-measured  distance  to  seaw  ater  and  the 
inslrument-detemiined  distance.  To  further  assess  this 
variation,  we  made  measurements  with  the  instrument 
along  a  1 .3-km-long  line  established  (not  forthis  study) 
across  first-year  sea  ice  of  varying  thickness  and  mor¬ 
phology  (Fig.  2).  Drill  hole  stations,  for  measuring  ice 
thickness  along  this  line,  were  spaced  5  m  apart.  EM3 1  - 
PM  instrument  soundings  were  made  at  206  of  the 


0  1  2 
H.  Tape-Measured  Distance  To  Seawater  (m) 

Fiifure  I .  EM31-PM  instrumcm-detcymincd  versus  tape-mea¬ 
sured  instrument  distanee  to  2.5-S/m-eonduciivity  seawater. 
\'ertieal  co-planar  coil  orientation  used. 


Distance  (m) 

Figure  2.  lee  relief  aloiif;  a  I  ..Fkm-lonf>  siin’ey  line  across  first-year 
sea  ire.  EM 3 1 -PM  instrument  ice  tliickne.ss  .soundinf’s  were  made 
alony  this  line 


stations  along  the  line  before  drilling  (F-ig.  3),  and  taped 
snow  plus  iee  thickness  measurements  were  made.  All 
EM3 1  -PM  soun  hng  measurements  were  made  with  the 
instrument  resting  on  the  snow  surface  as  shown  in 
Figure  4. 

A  plot  of  the  EM3 1  -PM  instrument-detennined  dis¬ 
tances  to  seawater  versus  the  drill-hole-measured  dis¬ 
tances  is  shown  in  Figure  5.  The  data  fall  into  two 
regions,  one  up  to  about  3.5  rn,  in  which  the  instrument- 
detennined  thicknesses  track  the  measured  distance  to 
seawater  reasonably  well,  and  a  second,  in  which  ex¬ 
tremely  poor  correlation  exists  for  o, stances  to  seawater 
of  over  3.5  m.  The  regression  curves  through  the  data 
representing  these  two  regimes  are  based  on  a  .some¬ 
what  arbitrary  3.5-m  break  point. 

The  w  interof  l'489-‘^0  produced  unusually  thick  -.ea 
ice.  Undeformed  first-year  sea  ice  with  0.05  m  or  less  of 
snow  cover  was  typically  2. 1  ±  0. 1  m  thick.  In  the  lee  of 
pressure  ridges,  snow  drift  depths  in  excess  of  1  m  we  c 
occasionally  encountered  on  the  level  sea  ice.  There¬ 
fore,  the  stand-off  distance  between  the  seawater  and 
the  EM3 1  -PM  on  undeformed  sea  ice  with  a  snow  cover 
could  reach  3  to  3.5  m.  The  regression  line  through  the 
0-  to  3.5-m-thick-ice  data  set  does  not  pass  through 
zero.  One  may  assume  that  this  is  caused  by  the  some¬ 
what  arbitrary  3.5-m  upper  bound  selected  for  the  data 
set.  as  well  as  the  paucity  of  data  at  the  higher  and  lower 


distances  or  ice  thicknes.s.  The  data  in  Figure  5  suggest 
that  EM3 1  -PM  sounding  can  provide  a  good  estimate  of 
snow  and  ice  thickness  from  about  0.7  to  3  m,  but  not 
with  an  accuracy  of  ±5%  of  the  directly  measured  value. 

Drill-hole-measured  snow  plus  ice  thicknesse.  over 
3.5  m  were  obtained  in  areas  of  deformed  ice.  The  poor 
agreement  between  the  EM3I-PM  instrument-deter¬ 
mined  distance  to  the  seawater  and  the  drill-hole-mea¬ 
sured  distance  in  these  areas  is  likely  attributable  to  the 
highly  variable  ice/water  interface  relief  in  the  area  of 
deformed  ice  and  pressure  ridges,  and  to  the  seawater- 
fdled  voids  in  the  ice  rubble.  These  voids  and  diffused 
interfaces  create  conductive  inhomogeneities  that  give 
rise  to  an  EMI  response  that  is  currently  not  inter¬ 
pretable.  Similar  unreliable  results  were  noted  by 
Hoekstra  (1979),  using  an  EM31  instrument,  and  by 
Kovacs  et  al.  ( 1 987b)  and  Kuvacs  and  Holladay  ( 1 989). 
evaluating  an  airborne  electromagnetic  induction  de¬ 
vice  for  sounding  sea  ice  thickness. 

After  reviewing  the  above  results  in  the  field,  and 
given  the  fact  that  there  was  no  thick  multi-year  sea  ice 
in  our  study  area  on  which  to  evaluate  the  EM3i-PM 
in.strument.  we  decided  to  replicate  thicker  ice  by  sim¬ 
ply  elevating  the  instrument  above  uniform  2.04-m- 
thick  sea  ice.  A  wooden  stepladder  was  used  to  elevate 
the  instrument  in  increments  up  to  about  2.8  m  above  the 
ice  surface  or  4.8  m  above  the  seawater. 


Figure  3.  Electric  drill  and  gu.v  engine  powered  units  used  to  drill  5-cm-diameter  ice 
thickness  holes.  The  EM31-PM  ele  tromai’netic  induction  soundint’  instrument  is 
restinf’  on  the  surface  in  the  foreiiround. 


D,  EM31 -Determined  Distance  to  Seawater  (m) 


li^ure  4 .  Typical  EM3 1  -PM  instrument  field  use  with  the  instrument  set  on  the  suifac 
wing  a  snow  plus  ice  thickness  measurement. 


Figure  5.  EM31-PM  instrument-determined  versus  tape-measi  '  "d 
distance  to  seawater  at  various  stations  along  a  1 .3-km-long  track 
across  level  and  ridged  first-year  sea  ice 


On-ice  calibration  of  the  EM3 1-PM  system  requires 
the  instrument  to  be  positioned  at  two  different  eleva¬ 
tions  ( H ]  and  hh )  above  the  seawater.  The.se  elevations, 
along  with  the  conductivity  of  the  seawater,  were  manu¬ 
ally  entered,  via  thumbwheel  dials,  into  the  PM.  The 
calibration  procedure  was  done  at  several  instrument 
elevations  to  detemiine  if  the  height  at  w  hich  the  in.stru- 
ment  was  calibrated  affected  the  EM,”^  1  -PM  instrument- 
detennined  ice  thickness  or  the  distance  to  the  seawater. 
In  addition,  for  one  calibration  a  value  of  3.5  S/m  was 
input  for  the  conducti\  ity  of  the  seawater  under  the  ice 
versus  the  true  value  of  2.5  S/m  as  detemiined  by  use  of 
an  in-situ  conductivity  probe. 

The  instrument  calibrations  were  made  at  H\  and  Ih 
distances  of  2.04  and  2.89  m.  2.04  and  3.45  m.  and  3.16 
and  4.02  m  using  the  correct  seawater  conductivity  of 
2.5  S/m.  For  the  case  w  here  3.5  S/m  was  used  for  the 
seaw  aterconducti  vity.  the  // 1  and  Hi  calibration  heights 
were  2.04  and  3.45  m  respectively. 

After  the  instrument  w  as  calibrated,  the  unit  was  set 
on  the  ice  surface  and  a  distance-to-seawater  measure¬ 
ment  made.  This  was  follow  ed  by  setting  the  instrument 
on  successise  steps  of  the  ladder  and  repeating  the 
distance  measurement.  The  resulting  data  are  listed  in 
Table  I  and  graphically  shown  in  Figure  6. 

The  plots  in  Figure  6  indicate  that  different  H\  and 
fh  instrument  calibration  heights  have  a  variable  effect 


on  the  instrument-determined  distance  to  the  seawater. 
While  the  spread  in  the  EM3 1  -PM  distance  determina¬ 
tions  for  the  various//;  and //2  calibration  heights  i.---  an 
the  order  of  ±5%,  in  most  all  cases  the  resulting  dis¬ 
tances  are  greater  than  the  tape-measured  distance  to  the 
seawater,  particularly  at  distances  over  3.5  m.  For 
example,  when  me  true  seawater  conductivity  was  used 
and  the  instrument  was  elevated  4.59  m  above  the 
seawater,  the  instrument  gave  distances  of  5.38,  4.88 
and  5. 13  m  for  test  runs  A,  B  and  C  respectively  (Table 
1).  The  average  of  these  reading  is  5.13  m  or  0.54  m 
greater  than  the  tape-measured  distance  of  4.59  m,  a 
13%  difference.  The  spread  not  only  becomes  larger 
with  increasing  distance  from  the  seawater  but  none  of 
the  regression  lines  passing  through  the  data  sets  inter¬ 
cept  zero  as  shown  in  Figure  7.  It  would  appear  that  the 
look-up  table  and  interpretation  algorithm  are  not  prop¬ 
erly  analyzing  the  received  electromagnetic  response 
from  the  seawater. 

It  is  interesting  to  compare  the  test  B  and  D  data  in 
Table  1.  Both  data  sets  were  collected  with  the  instru¬ 
ment  calibrated  at  the  same  H\  and  Hi  elevations,  but 
different  values  were  used  during  the  calibration  proce¬ 
dure  for  the  seawater  conductivity;  2,5  S/m  for  test  B 
versus  3.5  S/m  for  test  D.  A  plot  of  the  data  in  Figure  8 
shows  a  virtual  one-to-one  agreement  between  the  test 
results.  This  implies  that  for  the  test  stand-off  distances 


Table  1.  Laddertest  data  giving  the  tape-measured  distance  from  the  EM3 1-PM  to  the  ice  surface  and  ice/water 
interface  and  the  instrument-determined  distance  to  the  seawater  fur  different  calibration  heights  and 
seawater  conductivities, 
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Figure  6.  EM31-PM  instrument-determined  versus  tape-measured 
distance  to  seawater.  Data  were  obtained  by  elevating  the  instrument 
in  increments  on  a  stepladder.  Each  test  represents  a  different  instru¬ 
ment  calibration  height. 


H,  Tape- Measured  Distance  to  Seawater  (m) 


Figure  7.  A  replot  of  the  data  in  Figure  6  to  show  that  regre.'tsion  lines 
passing  through  each  test  data  set  do  not  intercept  zero. 
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Conductivity  Seawater  (m) 


Figure  8.  EM3I-PM  instrument-determined  distance  to  seawater  when 
3.5  S/m  was  set  into  the  instrument  for  the  water  conductivity  versus  the 
EM3J  -PM  instrument-determined  distance  to  seawater  when  the  correct 
conductivity  of  2.5  S/m  ua5  used  for  the  water. 


the  instrument-detennined  distance  to  seawater  was  not 
affected  by  using  the  above  seawater  conductivities  in 
the  calibration  procedure.  This  finding  will  not  hold  if 
the  instrument  is  brought  progressively  closer  to  the 
seawater.  When  this  occurs  the  response  from  the  sea¬ 
water  becomes  exponentially  larger  and  more  depen¬ 
dent  upon  the  water’s  conductivity.  It  is  assumed  that 
the  processor  module’s  ice  thickness  interpretation  al¬ 
gorithm  takes  this  effect  into  account. 

EM31  CONDUCTIVITY  READING  VERSUS 
SEA  ICE  THICKNESS 

As  previously  stated  the  electromagnetic  response 
measured  by  an  EM31  instrument  is  strongly  related  to 
the  conductivity  of  the  .seawater  and  the  stand-off  dis¬ 
tance.  The  conductivity  of  homogeneous  sea  ice  is 
about  two  orders  of  magnitude  less  than  that  of  Arctic 
Ocean  seawater,  and  should  have  a  minimal  effect  on 
the  electromagnetic  field  response  measured  by  the 
instrument.  Therefore,  we  decided  to  repeat  the  ladder 


measurements  but  only  record  the  EM3 1  ’s  conductivity 
reading  versus  instrument  distance  above  the  seawater. 
The  results  are  graphically  shown  in  Figure  9.  The  data 
clearly  show  that  a  standard  off-the-shelf  EM3 1  can 
provide  conductivity  measurements,  in  millisiemens 
per  meter,  that  are  directly  related  to  the  instrument 
height  above  seawater  or  to  ice  thickness. 

To  further  assess  this,  the  Geonics  PCLOOP  multi¬ 
layer  computer  code  was  used.  This  program  allowed 
calculation  of  the  response  that  an  EM3 1  instrument 
should  measure  versus  height  above  seawater  or  when 
the  instrument  is  resting  on  sea  ice.  The  program  uses 
operator  input  values  for  the  seawater  and  sea  ice 
conductivity,  as  well  as  the  distance  of  the  EM3 1  above 
the  seawater,  to  calculate  the  apparent  conductivity  that 
the  instrument  should  measure  under  these  conditions. 

In  this  analysis  the  seawater  conductivity  was  set  at 
2.5  S/m  and  a  2.04-m-thick  sea  ice  cover  with  a  bulk 
conductivity  of  1 0  mS/m  was  used.  A  comparison  of  the 
instrument-determined  field  reading  versus  the  theo¬ 
retical  apparent  conductivity  values  can  be  seen  in 
Table  2  and  in  the  plot  of  these  values  versus  the 
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H,  Tape-Measured  Distance  to  Seawater  (m) 

Figure  10.  EM3 1 -measured  and  computer-code  determined  conduc¬ 
tivity’  versus  tape-measured  instrument  distance  to  2.5-Slm-conduc- 
tivity  seawater.  Also  shown  is  the  difference  hetn’een  the  measured 
and  calculated  conductivity  values  versus  distance. 


tion  here  is  that  either  the  EM3 1  conductivity  determi¬ 
nations  are  wrong  or  the  computer-program  calculated 
instrument  distances  to  the  seawater  are  in  error.  From 
this  analysis  we  found  that  the  computer  code  consis¬ 
tently  calculated  lower  instrument-seawater  stand-off 


distances  than  determined  by  tape  measurement.  In 
addition,  subtracting  column  1  from  column  4  in  Table 
3  gives  the  results  in  column  5,  which  indicate  a  decreas¬ 
ing  ice  thickness  versus  increasing  instrument  height 
above  the  seawater.  A  plot  of  these  findings  is  shown  in 


Table  3,  Computer-code-calculated  distance  to  seawater  and  sea  ice  thickness  using  the  EM31  field- 
determined  apparent  conductivities  versus  instrument  elevation  above  the  seawater. 
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H,  Tape-Measured  Distance  to  Seawater  (m) 

Figure  II.  Ice  thickness  and  EM3 1  instrument  stand-off  distance 
to  seawater  determined  from  EM31  measurements  versus  tape- 
measured  distance  to  the  seawater. 


Figure  II.  Clearly,  there  is  a  problem.  The  problem 
turns  out  to  be  a  relatively  simple  one,  but  an  extremely 
important  one  that  must  be  appreciated  by  anyone 
choosing  to  use  an  EM31  type  instrument  for  sea  ice 
thickness  determinations. 

Resting  the  instrument  on  the  ladder  steps  caused 


two  things  to  happen.  First,  the  coil  support  booms, 
extending  out  each  side  of  the  electronic  module,  butted 
up  against  the  ladder  sides.  This  prevented  the  elec¬ 
tronic  module  from  resting  flat  on  the  steps  (see  Fig.  1 2). 
Second,  above  about  the  fourth  step,  the  tilt  was  exac¬ 
erbated  by  the  narrowing  of  the  sides  of  the  ladder  and 


Figure  12.  Position  ofEM3J  instrument  when  set  on  ladder  step. 
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by  the  operator,  who  stated  he  had  a  tendency  to  tilt  the 
unit  on  the  higher  steps  to  make  viewing  the  readout 
easier.  In  short,  he  did  not  wish  to  climb  any  higher  on 
the  relatively  unstable  ladder  than  was  apparently  nec¬ 
essary.  This  was  probably  a  good  idea,  considering  the 
poor  condition  of  the  stepladder  used. 

Tilting  of  the  EM31  ca  ised  the  orientation  of  the 
vertical  co-planar  coils  to  change.  Instead  of  being 
venical.  as  desired  and  assumed  in  the  computer  analy¬ 
sis,  the  coils  were  tilted  .shghtly  aw'ay  from  vertical. 

Fortunately,  the  computer  code  allows  detemiina- 
tion  of  the  theoretical  conductivity  for  both  a  vertical 
and  horizontal  co-planar  coil  position  for  a  given  set  of 
input  parameters,  such  as  instrument  height,  seawater 
conductivity  (2.5  S/m),  ice  thickness  (2.04  m),  ice 
conductivity  (10  mS;m),  etc.  Therefore,  the  computer 
code  was  used  to  recalculate  the  theoretical  conduc¬ 
tivity  for  the  vertical  and  horizontal  co-planar  coil 
positions  above  seawater  at  which  our  EM?1  field 
measurements  were  made.  These  calculations  are  given 
in  Table  4  (columns  2  and  3)  along  with  the  EM31 
reading  obtained  in  the  field  (column  4).  Next,  tilt  angle 
corrections  were  applied  to  the  calculated  values  in 
columns  2  and  3  until  an  angle  was  found  that  produced 
a  theoretical  conductivity  that  was  in  close  agreement 
with  the  field  measurement.  For  each  instrument  height, 
this  new  value  was  determined  by  multiplying  the 
number  in  column  3  by  the  cosine  of  some  coil  tilt  angle 
0  and  adding  this  value  to  the  value  obtained  by  multi¬ 


plying  the  number  in  column  2  by  the  sine  of  the  same 
angle. 

The  results  of  this  process  are  given  in  columns  5  and 
6  in  Table  4.  As  expected,  when  the  instrument  was 
resting  on  the  ice  surface  and  was  reasonably  level,  it 
gave  a  conductivity  reading  (column  4)  in  excellent 
agreement  with  the  theoretical  value  shown  in  column 
3.  However,  when  the  instrument  was  placed  on  the 
ladder  step  it  was  tilted  slightly  to  fit  securely  in  place 
(Fig.  12).  As  the  instrument  was  moved  onto  higher 
ladder  steps  the  tilt  angle  became  progressively  greater, 
as  shown  by  column  5  in  Table  4. 

Replotting  the  above  corrected  data  in  Figure  1 1 
would  cause  the  lines  to  rotate  counterclockwise  about 
the  2.04-01  intercept,  and  the  calculated  ice  thickness 
curve  would  became  horizontal,  as  it  should  be.  There¬ 
fore,  good  EM31  conductivity  measurements  should 
provide  for  a  determination  of  sea  ice  thickness  either 
for  when  the  instrument  is  resting  on  the  ice  surface  or 
when  it  is  elevated  above  the  ice  surface. 

Seawater  conductivity  will  be  lower  at  a  coastal  site 
where  a  major  river,  such  as  the  Mackenzie  River  in 
Canada,  continues  to  flow  into  the  ocean  during  the 
winter.  No  rivers  flow  into  the  Alaskan  Beaufort  Sea 
during  the  winter. 

As  previously  stated,  an  EM31  instrument  provides 
an  apparent  conductivity  reading  for  a  given  instrument 
height  above  seawater.  This  reading  is  unique  for  a 
specific  seawater  conductivity  and  should  allow  deter- 


Table  4.  Theoretical  versus  measured  EM31  conductivity  as  a  function  of  instrument  height  above  2.04-m- 
thick  sea  ice  and  instrument  tilt  angle  correction  required  to  produce  agreement  between  the  theoretical 
conductivity  for  the  vertical  co-planar  coils  on  the  EM31  and  the .  ctual  instrument  field  reading. 
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*  New  theoretical  vertical  coil  conductivity  value  =  Ver.  coil  conductivity  x  cos0  +  Hor.  coil  conductivity  x  sinS  . 
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Figure  13.  Theoretical  EM31  conductivity  reading  ver¬ 
sus  instrument  height above2.5-  and 3 .0-Slm  conductiv¬ 
ity  seawater. 


mination  of  instrument-seawater  stand-off  distance. 
During  the  winter,  the  seawater  under  the  Arctic  Ocean 
pack  ice  has  a  conductivity  of  2.5  ±  0.05  S/m.  An 
analysis  was  made  to  show  how  an  EM3rs  conductiv¬ 
ity  reading  should  vary  with  elevation  above  seawater 
of  2.5  S/m  and  an  extreme  seawater  conductivity  of  3.0 
S/m.  The  results  are  shown  in  Figure  13  for  vertical  co- 
planar  antenna  coils.  For  heights  above  about  3  m,  there 
is  no  appreciable  difference  in  instrument  reading  ver¬ 
sus  seawater  conductivity;  there  is,  however,  progres¬ 
sively  greater  difference  as  the  instrument  is  brought 
closer  to  the  water.  For  example,  when  the  instrument  is 
1  m  above  seawater  having  a  2.5-S/m  conductivity,  the 
EM31  reading  would  be  about  520  mS/m.  However, 
this  same  reading  over  3.0-S/m-conductivity  water  would 
occur  when  the  instrument  was  about  1.1m  above  the 
surface.  This  10%  error  for  the  extreme  seawater  con¬ 
ductivities  used  indicates  that  the  very  low  offshore 
Arctic  Ocean  conductivity  variations  will  not  signifi¬ 
cantly  affect  the  ice  thickness  determinations. 

Another  series  of  calculations  was  made  to  show 
how  an  EM3 1  ’s  conductivity  reading  should  vary  as  the 
instrument  is  moved  from  the  surface  of  1 .5-  and  2.5-m- 
thick  sea  ice  to  some  elevation  above  ice  floating  on 


seawater  having  a  conductivity  of  2.5  S/m.  In  this 
analysis,  bulk  sea  ice  conductivities  of  10  and  50  mS/m 
were  used.  The  latter  is  a  high  value  that  may  represent 
warm,  high-brine-volume  sea  ice.  The  former  value  is  a 
reasonable  value  for  arctic  winter  pack  ice.  The  results 
are  plotted  in  Figure  14.  As  may  be  inferred  from  these 
figures,  the  error,  caused  by  the  two  bulk  sea  ice 
conductivities  used,  in  the  distance  to  the  seawater,  as 
determined  from  the  instrument  conductivity  reading, 
is  on  the  order  of  5%.  This  indicates  that  the  response 
from  the  seawater  as  measured  by  the  instruments  will 
not  be  significantly  affected  by  sea  ice  conductivity 
variations. 

An  EM3 1  instrument  is  best  used  resting  on  the  ice 
surface  to  avoid  potential  measurement  error  associated 
with  tilting  of  the  antenna  coils.  This  would  be  particu¬ 
larly  desirable  on  thick  ice  to  maximize  the  received 
signal  and  to  allow  thick  ice  to  be  measured.  If  the 
instrument  must  be  carried,  then  some  provision  should 
be  made  to  show  the  operator  whether  or  not  the 
instrument  is  level. 

If  the  instrument  conductivity  reading  is  to  be  used  to 
infer  ice  thickness,  then  a  simple  table  or  graph  affixed 
to  the  instrument  cover  could  be  used.  Such  a  graph  is 
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H,  Instrument  Height  Above  Seawater  (m) 
a.  1 .5 -m-thick  ice. 


h.  2.5-m-(hick  ice. 


Figure  14.  Theoretical  EM3J  conductivity  reading  versus  in¬ 
strument  height  above  2.5-Slm-conductivity  seawater  covered 
with  ice  of  10-  and  50-mS/m  bulk  conductivity  (Gj). 
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shown  in  Figure  1 5.  This  figure  shows  how  the  conduc¬ 
tivity  reading  of  an  EM3 1  instrument  should  vary  with 
sea  ice  thickness  where  the  underlying  seawater  has  a 
conductivity  of  2.5  S/m.  Another  parameter  used  in 
construction  of  the  graph  is  a  variable  bulk  sea  ice 
conductivity.  These  bulk  conductivities  versus  ice 
thickness  were  obtained  from  recent  work  of  Kovacs 
( 1991 ).  Note  that  the  upper  curve  applies  to  an  instru¬ 
ment  resting  on  the  surface,  while  the  lower  curve 
represents  the  instrument  conductivity  readings  when 
no  ice  layer  exists 

The  sea  ice  conductivities  shown  in  Figure  1 5  repre¬ 
sent  reasonable  upper  limit  values  for  the  given  ice 
thickness.  Figure  15  is  instructive  because  it  indicates 
that  the  error  associated  with  ignoring  sea  ice  conduc¬ 
tivity  and  assuming  that  an  air  layer  exists  between  the 
instrument  and  the  seawater  will  have  a  minimal  effect 
on  the  estimated  ice  thickness.  This  is  especially  true  for 
ice  under  3  m  thick. 


GENERAL  COMMENTS 

Some  of  the  data  obtained  with  the  Geonics  EM31 
and  the  Flow  Research  ice  thickness  processor  module 
showed  excessive  variation  from  the  drill-hole-mea- 
sured  ice  thickness.  However,  measurements  on  first- 


Figure  15.  Theoretical  EM31  conductivity  reading  ver¬ 
sus  instrument  height  above  sea  water.  The  ice  is  resting 
on  water  of2.5-S/m  conductivity  and  ice  conductivity  is 
shown  to  vary  with  thickness. 


year  sea  ice  from  about  0.7  to  about  3.5  m  thick,  with  a 
snow  cover,  indicated  that  the  EM31-PM  instrument 
was  providing  reasonably  good  estimates,  even  if  they 
weren’t  accurate  enough  for  our  research.  The  unsatis¬ 
factory  ice  thickness  measurements  were  generally 
obtained  in  areas  of  defomied  ice  where  the  measured 
electromagnetic  response  was  adversely  affected  by 
conductive  inhomogeneities  associated  with  the  sub¬ 
merged  ice  block  structure. 

Further  testing  of  the  EM31-PM  instrument  on  a 
ladder  produced  ice  thickness  results  that  varied  with 
instrument  calibration  procedure.  This  should  not  oc¬ 
cur.  However,  other  inconsistencies  between  instru¬ 
ment-determined  sea  ice  thickness  and  direct  tape- 
measured  distances  may  have  resulted  from  tilting  of 
the  antenna  coils.  Given  this  mixed  review,  another 
field  test  of  the  module’s  measurement  performance  is 
needed. 

Use  of  the  EM31-PM  instmment  in  the  field  proved 
frustrating.  The  instrument  frequently  quit  at  tempera¬ 
tures  below  -15°C.  We  also  encountered  lithium  bat¬ 
tery  problems.  On  several  occasions  one  battery,  out  of 
the  set  of  ten  used  to  run  the  system,  would  drop  in 
voltage.  This  drop  would  shut  down  the  unit  until  the 
bad  battery  was  found  (back  in  “town”)  and  replaced. 
This  low  battery  problem  caused  measurement  delays 
and  may  also  have  caused  the  cold  weather  shutdown 
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problem  mentioned  above.  When  the  instrument  was 
picked  up  in  Kent,  Washington,  new  batteries  for  the 
instrument  were  also  obtained.  Therefore,  the  batteries 
should  have  been  good  if  they  were  recently  manufac¬ 
tured  ones. 

There  i s  good  reason  to  bel ieve  that  a  standard  EM3 1 
instrument  can  be  used  to  determine  sea  ice  thickness 
from  about  0.7  to  5  m,  to  within  an  accuracy  of  ±10%  of 
the  tape-measured  value.  For  this,  the  instrument’s 
conductivity  reading  and  a  table  or  graph  could  be  used; 
such  a  graph  is  presented  in  Figure  15.  It  may  well  be 
possible  to  achieve  ice  thickness  determinations  within 
5%  of  the  actual  thickness  when  on  ice  of  uniform 
thickness.  However,  further  field  tests  would  be  needed 
to  verify  this. 

In  the  absence  of  a  theoretical  analysis,  the  FM31 
may  be  calibrated  on  the  sea  ice  by  simply  recording  the 
instrument  conductivity  reading  made  at  a  number  of 
heights  above  the  ice/water  interface  and  constructing  a 
graph  similar  to  that  presented  in  Figure  4.  This  calibra¬ 
tion  procedure  needs  only  to  be  done  once  for  arctic  sea 
ice  sounding,  because  the  conductivity  determined  by 
the  FM3 1  is  based  on  the  use  of  the  measured  quadra¬ 
ture  pha.se  response,  the  zero  level  of  which  "is  designed 
to  stay  at  the  correct  value  over  the  life  of  the  instrument 
over  wide  temperature  excursions."*  However,  a  cali¬ 
bration  check  would  certainly  be  appropriate  at  a  loca¬ 
tion  where  under-ice  water  conductivity  may  not  be 
known. 

In  this  report  we  addressed  the  use  of  FMI  sounding 
to  measure  sea  ice  thickness  in  the  Arctic  Ocean.  The 
presentation  is  specific  to  areas  where  the  under-ice 
water  depth  is  in  excess  of  6  m.  In  shallower  waters,  the 
response  from  the  third  layer,  the  seabed,  needs  to  be 
addressed.  This  is  beyond  our  scope  here. 

It  should  be  kept  in  mind  that  FMI  sounding  does  not 
give  a  prist  measurement,  but  provides  an  "average" 
ice  thickness  for  an  area  having  a  diameter  of  about 
three  times  the  1. eight  of  the  instrument  above  the 
seawater.  Therefore,  FMI  sounding  of  multi-year  sea 
ice  thickness  may  result  in  significant  errors  where  the 
bottom  of  the  ice  has  rough,  short-period  undulations. 
This  relief  is  not  uncommon  and  must  be  taken  under 
consideration  when  multi-year  .sea  ice  soundings  are 
made  using  FMI. 


*  Personal  communication  with  J.D.  McNeill,  Gconics, 
Ltd..  1990. 
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